Full support of the Kasteleyn operator 
associated with a bipartite toroidal graph 
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Abstract. A perfect matching in a bipartite graph embedded on a torus defines a height function on the graph's 
faces and an associated height change vector in 1? . These matchings are enumerated by a combination of four 
evaluations of a bivariate Laurent polynomial, called Kasteleyn operator, whose coefficient of bidegree is, 
up to the sign, the number of perfect matchings with height change Therefore the Newton polygon of 

the Kasteleyn operator is the convex hull of the height change vectors. In this article, we prove that any point 
with integer coordinates in that polygon is realized by a perfect matching. 

1 Introduction 

Consider the set of doubly periodic lozenge tilings of the plane, for a given period C — (u|v)_BZ 2 , where B is 
an integer 2x2 matrix with nonzero determinant and u, v are the (column) vectors depicted in Figure la. For 
instance, the lozenge tiling shown in Figure lb has the period Z(3u + 4v, 5u + 4v). 




Figure 1 



The natural identification of lozenge tilings with stacks of cubes induces a height function that labels the 
points of Co := Z(u, v). Being the considered tilings invariant through translations by vectors of the sublattice 
C, such a translation produces a constant increment in the height function. For example, in Figure lb, the 
vector 3u + 4v reduces the height by 2 and the vector — 5u + 4v keeps it unchanged, so we can associate the 
pair (—2, 0) to this tiling. Not all the tilings with this period produce the same height change. For instance, 
the three constant tilings consisting of lozenges of the same type {a , w, and 0, respectively) give the following 
height changes: (1,9), (—11,-3), and (10,-6). 

For a general period C C £ (and a fixed basis), the possible height change vectors can be determined. Let 
h x ,h y , h z be the vectors associated to the constant tilings. If, for a ^-periodic tiling, the amount of lozenges of 
each type in a fundamental cell of the period is (x, y, z), its height change is the following convex combination: 
V °oi(C) ( x h-x + yh y + zh z ) — | dct(B) \ i x hx + yh y + zh z ), which is a point in the triangle determined by h x , h y , 
and h z - Moreover, its difference to any of these vertices lies in 3Z 2 , for the height difference between two points 
of Cq is determined modulo 3, independently of the tiling or the period. 

Conversely, for any point h in the triangle with coordinates in h x + 3Z 2 , there exists an ^-periodic tiling 
whose height change vector is h. The goal of this article is to give a proof of this fact in a more general setting, 
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which considers perfect matchings in a general bipartite graph embedded on a torus instead of periodic lozenge 
tilings. 

A slight modification of the construction explained above is a particular case of the height function defined 
by Kenyon et al. [2006] in the general case (see also [Thurston, 1990]). Summing up, a pair of perfect match- 
ings (loq,uj) defines a set of circuits (called transition graph) on the torus whose homology type h UQ (ui) £ 1? 
determines the height change vector. Details follow in Subsection 2.1. 

For a given bipartite toroidal graph and a fixed base perfect matching wqi the associated Newton polygon 
is the convex hull of the set of homology types h Uo (uS) of the transition graphs formed by u>q. It receives this 
name because it coincides with (a translate of) the Newton polygon of a bivariate Laurent polynomial P(w, z), 
namely, the determinant of a weighted Kasteleyn-Percus matrix (see [Kenyon et al., 2006]). The absolute value 
of a coefficient of P is the number of perfect matchings with a corresponding height change. In Section 3 we 
prove that any point with integer coordinates in the Newton polygon is realized by some perfect matching, so 
that the support (i.e. the set of occurring monomials) of the Kasteleyn operator is maximal. 

2 Preparations 

In this section we collect some basic facts that will be needed in the proof of our theorem. Firstly, we explain the 
formalization of the height functions described in the introduction, following [Kenyon, 2009]. Then, we discuss 
the homology of a set of knots in the torus. Finally, we relate the existence of circuits in circulant digraphs to 
the visibility of lattice points. 

2.1 Height functions 

In order to fix the class of graphs whose perfect matchings are assigned height functions, let us recall some facts 
about topological graphs. Wc refer to [Gross and Tucker, 1987; Mohar, 1988] for a detailed treatment. A graph 
G = (V, E) (possibly with loops or multiple edges) is endowed with a natural topology. An embedding of a graph 
G on a surface S is a mapping i : G c — > S such that the restriction i : G c — > is a homeomorphism (with 
i{G) endowed with the subspace topology). We usually identify a graph and its image through an embedding, 
so that G C S. Note that, as we deal with infinite graphs, the usual definition of embedding as an injective and 
continuous mapping is not equivalent. An embedding is cellular if the complement of the graph in the surface 
is homeomorphic to a disjoint union of open discs. 

We use the term periodic graph for a cellular embedding of a graph on the plane that is periodic through 
integer-valued vectors, i.e. the translation by any vector in Z 2 is a graph automorphism. The existence of a 
cellular embedding on the plane implies that the graph is locally finite. We will consider bipartite periodic 
graphs, but will not require them to be connected. 

The projection of a periodic graph G on the torus T 2 = R 2 /Z 2 is the embedding of a graph G. This graph 
is finite, for it is locally finite in a compact surface. Note that the embedding need not be cellular, as is shown 
in Figure 2. 

We use the notation A4(G) for the set of perfect matchings or 1-f actors of a graph G, i.e. 

M{G) := {u C E | Vu £ V 3!e € u : v e e}. 

In the following, we use the terms matching and perfect matching indistinctly, for we do not consider incomplete 
matchings. We are interested in those matchings on G which are compatible with the projection on the torus; 
or equivalently, matchings on G. 

For a bipartite periodic graph, we set a fixed orientation (from the white to the black vertex) at every edge. 
Given two matchings u, oj' € Ai(G), the transition graph lj — ui' is composed by the oriented edges of u> and the 
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Figure 2 



reversed edges of lj' . The connected components of a transition graph are directed circuits [transition cycles, 
according to Kasteleyn [1963]) and pairs of vertices bidirectionally linked, which we discard. We assume that 
the matching set is nonempty and fix a base matching uiq E M(G). Let u € M.(G) and let us define, by means 
of the transition graph u — lo , a height function on the faces of G, i.e. the vertices of the geometric dual G* . 
Note that G* is connected, although it need not be locally finite (and therefore, it is not necessarily a cellular 
embedding on the plane). Firstly, we choose a base face and assign the height to it. If there is an edge between 
two faces which does not occur in the transition graph, both faces get the same height. If the dual edge that 
goes from a face F\ to a face F 2 is crossed from left to right by the transition graph, we set ft,(i*2) = h(F\) + 1. 




Figure 3 



This process consistently defines an integer- valued function on the faces of G (see an example in Figure 3) . 
To see this, consider the three free Abelian groups generated by the vertex set of G, the set of edges endowed 
with the orientation defined above, and the set of (oriented according to the embedding) faces, respectively. 
The sketch below represents the border homomorphisms <9; and their transposed: the coborder homomorphisms 
Si. 

5 61 



Z(V(G)) Z(E(G)) Z(F(G)) 

di d 2 



Writing d* and 5* for the border and coborder homomorphism (respectively) associated with the dual graph, 
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we have 81 = 5%, <9 2 — Sq. An element t* £ Z(E(G*)) defines the homomorphism 



t* = (r, •) : Z(E(G)) — >■ Z 

w 1 y ^r(e)w(e), 

eG-B 

and an analogue construction can be done for elements of Z(F(G*)). With this notation, if F\,F2 £ F(G) and 
t* is a walk from F\ to F 2 in the dual graph, we have h(F 2 ) — h(F 1 ) = (t,lo — loq). Now, if and t 2 are 
walks in G* with the same endpoints, we have to see that (t\,lo — loq) — (t 2 ,w — loq). Note that r* — r| is a 
border (there is p* £ Z(F(G*)) such that 8*(p*) = t{ - r|) and lo - lo is a cycle (<9(w — w ) = 0). Therefore 
(n -r 2 ,w- wo) = (S(p),lo - w ) = (8(uj - lo ),p) = 0. 

On the other hand, if F and F' are faces of G and u <E Z 2 , we have + u) - h(F) = ^(i^' + u) - /i(-F')- 
The height function is determined, therefore, by its values on a system of face representatives modulo Z 2 and 
the pair of height increments h^^Lo) :— (h(F + (1, 0)) — h(F),h(F + (0, 1)) — h(F)). As before, it can be shown 
that the height change vector is determined by the homology type of to — ui in the torus. It is easy to check 
that, if the homology type of u) — ujq is h U(i (<jj) = (a,b) £ Z 2 , we have /i Wo (w) — (—6, a). For instance, in the 
example of Figure 3, we have /^(w) = (1, —1) and h Uo (uj) = (—1, —1). 

2.2 Torus knots 

As we have seen, the height change of a matching (with respect to a base matching) can be identified with the 
homology type of a set of disjoint oriented copies of S 1 on the torus. Let us collect some results on torus knots 
(whose proofs can be found in [Rolfscn, 1976]) for later use. Note that, unlike usually, we consider directed 
knots. 

Lemma 1 Let u = (111,112) £ Z 2 . There exists a torus knot c with homology type u if and only if u is a visible 
lattice point (i.e. u = or gcd(ui,u 2 ) = 1). 

Let / £ Aut(T 2 ) be a self-homeomorphism on the torus. The functor 7Ti associates an automorphism of the 
fundamental group Z 2 with it: 

Aut(T 2 )^GL(2,Z). 

This mapping is a group epimorphism, in particular, for every group automorphism P on Z 2 , there is a self- 
homeomorphism on the torus whose effect on the homology types of the torus cycles is determined by P. 
Moreover: 

Lemma 2 Let c\ , c 2 be two torus knots with nontrivial homology. Then there is a self-homeomorphism on the 
torus which maps c\ into c 2 . 

Lemma 3 Let c\ , c 2 be two disjoint torus knots with nontrivial homology types. Then both homology types 
coincide or are opposite. 




1 1 
1 i 



Proof. By Lemma 2, we can assume that the homology type of c\ is (1,0). Therefore c 2 is a knot on the 
cylinder T 2 \ci, so its homology type must be ±(1, 0). ■ 
The following result is a direct consequence. 

Lemma 4 Let G be a bipartite periodic graph and wo, w £ A4(G), such that u = (11%, u 2 ) := h Un (uj) £ Z 2 \{0}. 
Set d :— gcd(wi, u 2 ). Then the transition graph consists of circuits with zero homology, P circuits with homology 
type ^u, and N circuits with homology type -jfU, with P — N = d. In particular, if the transition graph consists 
of a single circuit, u is a visible lattice point. 

Note that if a transition graph lo — lo q consists of several circuits, the removal of some of them leads to 
another transition graph lo' — loq. Using this remark, we get: 
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Lemma 5 Let G be a bipartite periodic graph anduiQ,ui G Ai(G), such that u = (ui,tt2) := /i Wo (u;) G Z 2 \{0}. 
5et d := gcd(tti, ^2)- Then there is a matching u 1 € M.(G) such that the transition graph u/ — uq consists of 
one single circuit with homology type ^u. 

2.3 Circulant digraphs 

In the proof of the main result of this article we use a certain class of toroidal graphs whose matchings can 
be identified with sets of disjoint circuits in directed circulant graphs. Let us recall the definition and some 
properties of these objects (see the survey by Bermond et al. [1995]). 

For integers n, j\, . . . ,j r such that n > 1, the circulant graph C(n;ji, . . - ,jt) is the Cayley graph with vertex 
group Z/nZ and edges {{x, x + ji} | 1 < i < t}, where the jumps ji are considered as residue classes modulo n. 




Figure 4: C(8;l,3) 



Note that we are including in the definition circulants with repeated jumps (leading to multiple edges) and 
zero jumps (leading to loops). We denote by C(n\j\, . . . , jt) the corresponding directed graph. For instance, 
the digraph C(n; 1, 1) consists of a directed cycle of length n with duplicated arcs. 

Let us consider a circulant digraph with two jumps: C(n;a, b). If a walk consists of u arcs of the type 
(x, x + a) and v arcs of the other type, we say that its abelianized is (it, v). If x G Z/roZ is the initial vertex of 
such a walk, its terminal vertex is x + au + bv. We can label, as is done in [Wong and Coppersmith, 1974], a 
point (u, v) G N 2 by au + bv: the terminal vertex of a walk with abelianized {u,v) and initial vertex 0. 

Note that the set of abelianized closed walks in a connected graph C(n;a, b) is the intersection of a 2-rank 
integer lattice with volume n (which we call circuit lattice) with N 2 . For instance, in the case C(n;l,6), the 
circuit lattice is Z((n, 0), (—6, 1)). 

We use the term lattice path for a finite list of points in N 2 such that the difference between two consecutive 
points is (1,0) or (0, 1). The labelling considered above allows the identification of walks in the digraph with 
lattice paths. When there are no two distinct points in a lattice path which are congruent modulo the circuit 
lattice, the associated walk in the digraph is a path (i.e. no vertex is visited twice). When the only pair of 
congruent points modulo the circuit lattice consists of the endpoints of the lattice path, the associated walk is 
a circuit. 

Lemma 6 Let A C Z 2 be an 2-rank integer lattice and v G A n N 2 . There is a lattice path from the origin to 
v such that no difference between two distinct path nodes ( except v — ) lies in A if and only if v is visible in 
A (i.e, the segment joining and v does not contain any other point in A) and ||v||i < vol(A). 

Proof. The condition ||v||i < vol(A) is necessary: the group Z 2 /A has vol(A) elements, a lattice path with 
length bigger than vol(A) visits at least vol(A) + 2 points, and therefore, at least two pairs of them are congruent 
modulo A. 

If v = (1)1,1)2) is not visible in A, the required lattice path cannot exist either. In that case, there is an 
integer d > 2 such that G A. Consider a lattice path p = (0 = po, ■ ■ ■ >P||v||i = v ) an d its subpaths p(i, d) 
with origin at pi and length h |jv|j 1, for i = 0, . . . , ^p||v||i. Let rij be the number of steps of type — > that p(i, d) 
consists of. If, for some index i, we have = v\/d, the difference of the endpoints of p(i, d) is ^ v G A. In other 
case, since Uj+i — n\ G {1, 0, —1}, we have: rii < v\/d, for every index i; or m > v\/d, for every index i. Then, 

d-i 

It can be interesting to compare this implication with the so-called Universal Chord Theorem (see [Rolfscn, 
1976, p. 15]). For the other, let v G A n N 2 be a point visible in A such that ||v||i < vol(A). Draw the segment 
s that joins the origin with v and consider the diagonals rf, := {(x,y) G M 2 | x + y = i}. We define the path 
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P = {POiPii ■ ■ ■ )P\\v\\i) by choosing pi as the element of di with integer coordinates that minimizes the distance 
to s n di (see Figure 5a). If the intersection is equally distant from two integer points, we choose the one with 
a bigger first coordinate. It is easy to see that p i+ i — p t € {(1, 0), (0, 1)} and p is a lattice path from to v. 




(a) (b) 

Figure 5 



Suppose that < i < j < ||v||i, j — i < ||v||i, and pj — pi E A. Note that pj — pi and v must be linearly 
independent, for \\pj — Pi\\i = j — i < ||v||i and v is visible in A. Consider the lattice path p' "parallel" to p 
with origin at pj — pi, i.e. p' k := p^ + pj — p il for k = 0, . . . , ||v||i. These paths collide, for p[ — pj. This implies 
that the image of A through the projection onto K((l, —1)) parallel to s has two points with distance smaller 
than (a pair of points distant exactly \/2 is not enough, as we have set up a rule for tie breaking). 

As v is visible in A, there exists w £ Z 2 such that A = Z(v,w). The projection described above is defined 
by the equation: 



so that 7r(A) = Z(-7r(w)) = Z(^y^(l, — 1)), and the distance between two projected lines is at least v °^\^ V%- 
Therefore the paths cannot collide and the proposed lattice path satisfies the requirements. ■ 

We close this section with a digression, obtaining, as a corollary of the previous result, the well-known 
characterization of circulant digraphs with two jumps which are Hamiltonian. Curran and Witte [1985] used 
similar arguments for the study of Hamiltonian paths. 

Consider a connected circulant digraph C(n;a, b), i.e. gcd(n, a, 6) = 1. A Hamiltonian circuit in this graph 
corresponds to a lattice path starting at the origin 0; ending at a certain v; and such that every element in 
TLjnL is the label of exactly one visited point, except 0, which labels both and v. In particular, ||v||i = n. As 
the volume of the circuit lattice of C(n; a, b) is n, Lemma 6 gives the following characterization of Hamiltonian 
circulant digraphs: 

Corollary 7 Let n be a positive integer and a,b € Z such that gcd(a, b, n) = 1. Let A := {(u, w) e Z 2 au + bv € 
nli}. The circulant digraph C(n; a, b) is Hamiltonian if and only if the intersection An{(x, y) £ N 2 | x + y = n} 
has a point visible in A. 

Note that the intersection of the circuit lattice with the diagonal {x + y = n} and the first quadrant is 
{(n — «§,*§) | i = 0, . . . , d}, where d := gcd(& — a, n). The point (n — «§,«§) is visible in A if and only if the 
label of (d — i, i) has order n/d in TLjnL. Therefore Corollary 7 is equivalent to the following characterization 
(see [Fiol and Yebra, 1988; Locke and Witte, 1999; Yang et al., 1997]), whose proof dates back to [Rankin, 
1948]. 

Proposition 8 Let n be a positive integer and a,b € Z such that gcd(a, b, n) = 1. The circulant digraph 
C(n;a,b) is Hamiltonian if and only if there exist nonnegative integers i,j such that: 

i + j = gcd(& — a, n) — gcd(ai + bj, n). 
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3 Main result 



As we have seen in Section 2.1, given a bipartite periodic graph G and a matching wo on its projection G, 
any matching ui G .M(G) defines a point h^uj) e Z 2 , related to its height function. Consider the associated 
Newton polygon N ulQ (G): the convex hull of those points. We prove that any point with integer coordinates in 
this polygon is h Uo (uj), for some uj G M(G). This problem can be reduced to visible points in a triangle with 
one vertex at the origin and the others at visible lattice points (named v and w) . An example of a graph whose 
Newton polygon attains this shape is depicted in Figure 6b. It is obtained by projecting the segments [0,v] 
and [0,w] onto the torus and substituting an u>o-edge for every crossing. 




(a) (b) 

Figure 6: v = (3,2), w = (1,4) 



For integers n, r such that < r < n, we consider the bipartite periodic graph B(n, r) with black vertices at 
points (i/n,j) and white vertices at points ((2i + l)/2n,j), for integers i,j. Each white vertex ((2i + l)/2n,j) 
has the three following neighbours (see Figure 7): (— , jj , , (^^i j + !)• The edges of the first type 

constitute a perfect matching ujq, which we fix as base matching. 




Figure 7: A projection of B(5, 2) 



Proposition 9 Let n,r £Z such that < r < n and gcd(r, n) = I, and consider the (closed) triangle T with 
vertices (0,0), (1,0), and (r,n). If u is a visible lattice point in T, there is a matching lu £ M{B(n,r)) such 
that h^ (ui) = u. 

Proof. Consider the projected graph B{n,r) and identify vertex (i/n,j) with ((2i + l)/2n, j), obtaining the 
circulant graph G(n;l,r). Matchings in B(n,r) correspond with sets of disjoint circuits in C(n;l, r), whose 
circuit lattice is A := Z((n, 0), (— r, 1)), and the following bijection relates the homology type of a matching to 
the abelianized of the corresponding set of circuits. 

t: TnZ 2 — > An {(x,y) € N 2 | x + y < n} 

(:) - (o t)(:) 

Let u = (u\, U2) 6 Z 2 be a visible point in TnZ 2 . By Lemma 6, there is a circuit in C(n; 1, r) whose abelianized 
is t(u). If w is the associated matching, the homology type of ui — uj is u. ■ 
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Theorem 10 Let G be a bipartite periodic graph and ujq € A4(G). We have: 



A^ (G)nz 2 = {h 

u 



,(w) | oj e M(G)}. 



Proof. We have to prove that the set on the left-hand side is contained in the other, which we denote by h(G). 
Let u S Nu (G) HZ 2 . It is easy to see that this element is in h(G) if it is a vertex of Nu (G) or u = = h^ (ojq). 

Note that, for any three matchings w^o^, w £ A4(G), /i Wl (u;) = h ull (u)2) + h u)2 (w). Therefore a change in the 
base matching is equivalent to a translation in the corresponding set h(G). Suppose that v,w g h(G). Using 
Lemma 4, any other point with integer coordinates in the segment they define lies also in h{G). 

By Caratheodory's theorem, u lies in the triangle defined by and two vertices of N UJo (G). Using the 
previous argument, we can assume w.l.o.g. that u is a visible point in the triangle defined by and two visible 
and linearly independent points v, w £ h(G). By Lemma 5, there exist uj±, UJ2 € M{G) with respective homology 
types v and w, and such that each of the transition graphs u)\ — ujq and UJ2 — consists of a single circuit. 

We know (see the discussion in Subsection 2.2) that there is a self-homeomorphism / on the torus whose 
associated automorphism on the first homology group satisfies: 



tti(/)(v|w) 



where < r < n := |det(v|w)| and gcd(r, n) = 1. Using Lemma 2, another torus self-homeomorphism 
transforms the transition cycle u)\ — ujo into the meridian [0, 1] x {0}. The set of height change vectors of the 
transformed graph G\ is 7r 1 (/)/i(G), so we need just prove that any visible point in the triangle defined by 
(0,0), (1,0), and (r,n) is in h{Gi). 

Now, we transform G\ into another graph H such that h(H) C h(G\). Consider the intersections of the 
meridian (image of w\ — u>o) an d the image of the circuit u>2 — ^o- Such an intersection cannot be limited to a 
single point. If this was the case, that point should be a vertex incident to two different edges of ojq, which would 




U) 



w 




u 




(b) 



Figure 8 



not be a perfect matching (see Figure 8a). Therefore each connected component of the intersection consists of 
an odd number of edges, which can be reduced to one, as depicted in Figure 8b. 

Suppose that in G\ there are two consecutive (as walking through the meridian) crossings between both con- 
sidered circuits with the same F-coordinate. We can remove (see Figure 9) all those crossings (which happen in 




Figure 9 



finite number, for the projected graph is finite), obtaining a graph H homeomorphic (as a planar embedding) to 
B(n, r). Indeed, after the removal process, there are exactly n crossings in H, and the edges of the second circuit 
define a permutation of them. As there is no self- intersection, this permutation must be of the type x i— > x + s. 
We have gcd(s, n) = 1 (s = r, indeed), for ll>2 — ojq consists of a single circuit. Therefore Proposition 9 applies. ■ 
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